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1.  INTRODUCTION 


The  geomagnetic  field  plays  an  important  role  in  a > * -^r  of 

high-altitude  nuclear  explosion  phenomena,  including  debris-air 
coupling,  the  guiding  of  beta  rays  and  energetic  ions  down  into  the  at- 
mosphere, and  the  formation  of  striations,  to  mention  a few.  For  first 
bursts,  and  sufficiently  late  after  any  burst,  this  field  will  be  the  am- 
bient geomagnetic  field.  Accordingly,  a model  of  the  ambient  geomag- 
netic field  is  needed  for  the  new  radar  and  optical  systems  code. 


The  requirements  of  such  a model,  apart  from  the  general  ones 
of  modularity  and  minimal  demands  on  computer  storage  and  running 
time,  are  that  it  provide  reasonably  accurate  values  of  the  vector  field 
components,  and  that  it  permit  the  efficient  tracing  of  field  lines. 

The  RANG  codes  used  an  eai'th- centered  dipole  approximation 
to  the  ambient  geomagnetic  field.  Such  a model  is  certainly  fast,  and 
permits  the  easy  tracing  of  field  lines,  but  its  predictions  are  of  low 
accuracy.  On  the  other  hand,  there  are  available  highly- accurate 
multipole  field  models  [SM-72e]  that  are  fast- running  except  for  their 
field- line  tracing  routines,  which  necessarily  integrate  numerically. 

Because  of  the  fact  that  the  systems  code  will  be  concerned 
with  only  a limited  ’ ittle  space  of  the  order  of  one  thousand  kilometers 
in  linear  dimension,  a compromise  solution  incorporating  the  best 
features  oi  both  kinds  of  model  becomes  possible;  it  has  been  explored 
and  is  tentatively  adopted.  This  model  uses  accurate  field  components 
obtained  from  the  multipole  model  for  some  point  in  the  middle  of  the 
battle  space  to  fit  a locally-best  geocentric  dipole  field.  This,  of 
course,  needs  to  be  done  only  once,  during  problem  setup.  The  dipole 
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model  is  then  used  for  subsequent  field  evaluations,  for  line  tracing, 
and  so  on.  Thus  both  speed  and  good  accuracy  are  obtained. 

In  the  following  sections  there  is  a description  of  a set  of  com- 
puter subroutines  that  have  been  written  to  implement  the  model. 
Listings,  cross-reference  lists  of  variables,  and  input/output  lists  are 
included  in  an  appendix,  along  with  test  problems  that  have  been  used 
in  model  evaluation. 


I 
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2.  SUBROUTINE  ONEMG5  AND  LINTRA 

Personnel  of  the  National  Aeronautics  and  Space  Administration 
hav»  developed  and  thoroughly  documented  [SM-72e]  a set  of  Fortran 
subroutines  providing  a multipole- expansion  model  of  the  geomagnetic 
field,  including  secular  changes  and  provisions  for  tracing  field  lines 
lO  intersects  at  specified  altitudes.  One  of  these  routines,  called 
ONEMG5,  * embodies  the  International  Geomagnetic  Reference  Field 
(IGRF  1965. 0),  and  it  has  been  adopted  here  as  the  "good”  magnetic 
field  model.  Another  routine  called  LINTRA  traces  geomagnetic  field 
lines  to  their  intersections  with  prespecified  altitudes;  it  has  been  used 
only  for  verification  of  the  simplified  dipole-field  line-tracing  routine. 

A simplified  flow  diagram  of  ONEMG5  is  shown  in  Fig.  1. 
Corresponding  details  for  LINTRA  have  not  been  supplied  here,  for  that 
routine  does  not  form  part  of  the  present  package,  but  was  only  used  in 
evaluation.  Moreover,  these  details  are  readily  available  in  SM-72e. 

One  note  of  warning  must  be  sounded  concerning  the  descrip- 
tion of  secular  changes  that  is  provided  in  CNEMG5.  This  description 
is  of  first  order  only,  and  is  based  on  a fairly  small  number  of  years  of 
good  data  near  the  epoch  1965. 0.  Consequently,  it  is  inadvisable  to 
input  a time  xnore  than  a few  years  away  from  the  data  range  of  the 
model. 


’•'Called  ONEMAG  in  £;M-72e. 


3.  SUBROUTINE  MAGFIT 

Given  a point  in  space  (normally  near  the  earth's  surface  and 
centrally  located  in  the  battle  space)  for  which  accurate  values  of  the 
geomagnetic  field  components  are  known,  subroutine  MAGFIT  calculates 
the  strength  and  orientation  of  an  earth-centered  magnetic  dipole  to  re- 
produce those  components.  The  routine  is  used  only  once,  during 
problem  setup,  and  the  dipole  properties  are  then  stored  and  used  later 
to  provide  field  component  values  at  other  points  within  the  limited 
battle  space. 

In  Fig.  2,  the  point  P at  geocentric  radial  distance  r,  north 
latitude  X (colatitude  6),  and  eart  longitude  (f>  is  the  reference  point  at 
Which  the  field  components  B^,  Bq,  and  B^  (in  the  same  coordinate 
.system  (r,  6,  ip)  are  known.  The  point  Q at  north  latitude  (colati- 
tude 6^)  and  east  longitude  <p^  on  the  surface  of  an  earth- centered 
sphere  passing  through  P is  the  direction  of  the  earth- centered 
dipole.  The  point  N is  the  north  geographic  pole.  The  arc  length  (or 
Central  angle)  between  Q and  P is  denoted  by  x- 

From  the  equations  for  a magnetic  dipole  field  we  have  the 

relations 

= , (1) 

®2  " ’ (2) 

Where  M measures  the  dipole  strength  and  Bg  is  the  angular  component 
in  the  direction  of  increasing  X.  Consequently,  from  simple  geometry 
there  follows  the  relaiicns 


and 
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B, 


B 


b: 


<p 

2 


= cos  a , 

(3) 

= Bg  sin  a , 

(4) 

- ^ 

(5) 

where  ot  is  the  angle  QPN. 

From  Eqs.  (1)  and  (2)  one  finds  the  formulas 
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X = tan"-  (2B2/Bj.)  . (7) 

From  Eqs.  (3)  and  (4)  it  follows  that 

a = tan"^  (B^Bq)  . (8) 

By  applying  the  cosine  law  of  sjAerical  trigonometry  to  the  spherical 
triangle  QPN,  one  obtains  the  relation 

cos  0 = cos  X cos  6 + sin  X sin  0 cos  a . (9) 

o 

Application  of  the  sine  law  leads  to  the  further  relation 

sin  (<p  - <f>^)  = sin  x sin  a/sin  0^  . (10) 

One  more  use  of  the  cosine  law  yields  the  equation 

cos  («p  - «p^)  = (cos  X - cos  0^  cos  0)/(sin  0^  sin  0)  , (11) 

useful  in  establishing  the  correct  quadrant. 
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Equations  (5)- (11)  constitute  the  working  equations  of  subroutine 
MAGFIT.  A Fortran  listing  of  the  routine  appears  in  the  appendix.  A 
simplified  flow  diagram  is  given  in  Fig.  3. 


Subroutine 

MAGFIT 


Set  up  input  parameters  for  call 
to  magnetic  multipole  field  model 


Calculate  magnetic  field 
components  at  point  P using  mag-  ■ 
netic  multipole  model  (IGRF  1965.0) 


I Compute  magnetic  dipole  moment 
I from  the  field  components 


RETURN 


Fig.  3.  Subroutine  MAGFIT  Flow  Diagram. 


Subroutine 

ONEMG5 


4.  SUBROUTINE  BFIELD 


For  any  point  P at  geocentric  radial  distance  r,  north  latitude 
X (colatitudc  6),  and  eas<'  longitude  '4>,  subroutine  BFIELD  calculates 
the  geomagnetic  field  strength  B,  the  dip  angle  I,  and  the  declination 
angle  D,  based  on  a locally-fitted  geocentric  magnetic  dipole  of  strength 
M oriented  in  the  direction  of  north  latitude  X^  (colatitude  6^)  and  east 
longitude  <p  . These  latter  three  quantities  must  have  been  found  previ- 
ously by  the  use  of  subroutines  MA^GFIT  and  ONEMG5  for  a reference 
point  within  a thousand  kilometers  or  so  of  point  P,  if  good  accuracy  is 
to  be  assured.  Figure  2 may  be  used  to  help  visualize  the  geoiinetrical 
relationships. 

The  equations  of  subroutine  BFIELD  are,  for  the  most  part, 
just  those  presented  above  in  Section  3,  but  solved  for  different  vari- 
ables. Thus,  the  angle  X between  the  dipole  moment  and  the  field  point 
is  obtained  from  the  equation 

cos  \ = cos  0 cos  0 + sin  0 sin  0 cos  {(p  - (p)  , (12) 

O 0 0 

which  follows  from  spherical  trigonometry.  The  total  field  strength  B 
is  obtained  by  use  of  Eqs.  (1)  and  (2)  from  the  relation 


= ^ [3  cos^x  + 1]^ 


The  dip  angle  I is  obtained  by  use  of  Eqs.  (1)  and  (2)  and  the 


definition 


sin  I = B /B 
r 


= 2 cos  x/  3 cos  X + 1 


f- 


I 


r( 


r’ 


5.  SUBROUTINE  CONJUG 

The  main  function  o!  subroutine  CONJUG  is  to  locate  the  lati- 
tudes and  longitudes  of  those  points  where  a given  geocentric  magnetic 
dipole  field  line  intersects  a prespecified  altitude.  There  are  generally 
two  such  points;  the  routine  will  locate  either,  depending  on  the  choice 
of  an  input  quantity.  CONJUG  also  computes  (1)  the  dimensionless 
field- line  distance  (in  units  of  the  equatorial  radius  to  the  dipole  field 
line)  between  two  specified  points  and  Pg  and  (2)  the  ratio  of  the 
equatorial  field  to  that  at  point  Pj^,  for  the  same  field  line. 

Suppose  the  orientation  of  the  geocentric  dipole  is  specified  by 
the  north  latitude  (colatitude  6^)  and  east  longitude  (p^.  Let  the 
field  line  be  specified  by  the  fact  that  it  passes  through  a point  P^  in 
space  at  altitude  h^,  north  latitude  (colatitude  0^),  and  east  longi- 
tude <py  Then  we  seek  the  north  latitude  Xg  (colatitude  0g)  and  east 
longitude  <P2  of  a point  P^  on  the  same  dipole  field  line  as  P^.  The 
geometry  of  the  situation  is  illustrated  in  Fig.  5. 

From  the  cosine  law  of  spherical  trigonometry  applied  to  the 
spherical  triangle  P^  P^  N,  we  obtain  the  result 

cos  Xj  = cos  0^  cos  0j  + sin  0^  sin  0^  cos  (ip^  - <P^)  • (20) 

The  sine  law  for  the  same  triangle  gives  the  result 

sin  0 = sin  0j  sin  (ipj  - <p^)/sin  Xj  » (21) 

and  another  application  of  the  cosine  lav/  gives  the  formula 

cos  0 = (cos  Gj  - cos  Xj^  cos  0^)/sin  Xj  sin  0^  , (22) 


so  there  is  no  ambiguity  as  to  the  quadrant  of  tj). 


f 
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The  equation  of  a dipole  field  line  has  the  form 

. 2 

r = sin  x f 


so  the  requirement  that  and  Pg  lie  on  the  same  dipole  field  line 
leads  to  the  result 


sin  Xg 


sinxj 


(24) 


where  R is  the  radius  of  the  earth  and  h^  is  the  prespecified  altitude 
of  point  Pg*  Note  that  there  are,  generally,  two  solutions  for  X2> 
since  if  X2  is  a solution,  so  is  ff  - X2* 

Now,  applying  the  sine  and  cosine  laws  to  spherical  triangle 

P«  P«  N leads  to  the  results 
o 2 

cos  6«  = cos  6 cos  Xo  + sin  6^  sin  Xq  cos  tp  (25) 

Z O 6 O 6 

and 

sin  ((pg  - <p^)  = sin  X2  sin  0/sin  ©g  , (26) 


whence  0g  (or  Xg)  and  ipg  can  be  obtained. 


The  absolute  value  of  the  dimensionless  field- line  distance  be- 
tween points  Pj  and  Pg  is 


S 


12 


/ 


ds 


where  the  element  of  arc  length  is  given  by 


^ = r^  sin  X (1  + 3 cos^x)^  . 
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(2 


'f) 


(28) 


After  substituting  Eq.  (28)  into  (27)  and  performing  the  integration,  we 
obtain 


^12  = f "i  Vi+'if  - \ ill  *4  * 


where 


7}^  = *73  cos  Xj 


(30a) 


t?2  = cos  Xg  • 


(30b) 


Equation  (29)  is  valid  provided  points  and  Pg  are  in  the  same  hemi- 
sphere. If  points  Pj  and  Pg  are  in  opposite  hemispheres,  then  we 
must  perform  the  integration  in  two  parts,  with  the  equator  being  the 
intermediate  point.  The  result  may  be  expressed  in  the  form 


5 = - AJUG  X Sr 

12  6 IE  ^5 


(31a) 


where 


®1E  “ "l  Vl*"!  * fl 

®2E  ' '>2  Vl  + ”2  (”2  + Vl  + ’>2) 


(31b) 


(31c) 


and  AJUG  is  a parameter  equal  to  (+1)  if  Points  P^  and  Pg  are  in  the 
same  hemisphere  and  equal  to  (-1)  if  Points  Pj  and  Pg  are  in  opposite 
hemispheres. 


The  equatorial  radius,  r^,  is  given  by 


To  = (M/Bf  , 


Fig.  3.  Subroutine  CONJUG  Flow  Diagram. 
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6.  SUBROUTINE  MAGDRV  AND  VERIFICATION  TESTS 
OF  THE  AMBIENT  GEOMAGNETIC  FIELD  MODEL 


To  permit  the  exercise  of  the  ambient  geomagnetic  field  model 
for  purposes  of  testing  and  validation,  a special  driver  routine  called 
MAGDRV  has  been  v;ritten.  The  required  input  consists  of  latitude, 
longitude,  and  altitude  coordinates  of  a set  of  reference  locations,  at 
each  of  which  the  vector  field  of  a geocentric  magnetic  dipole  is  fitted 
to  an  accurate  multipole  field,  for  a specified  year.  Further  input  con- 
sists in  a set  of  locations,  relative  to  each  reference  point,  for  which 
both  the  dipole  field  and  the  accurate  multipole  field  are  evaluated  and 
compared  for  relative  accuracy  of  the  total  field  strength.  Additional 
input  consists  in  sets  of  altitudes  for  the  calculation  of  field- line  inter- 
sects for  each  of  the  test  points,  together  with  flags  indicating  whether 
the  desired  intersection  is  in  the  same  or  opposite  magnetic  hemisphere. 
Additional  output  consists  in  the  inclination  and  declination  angles  for 
each  test  point,  according  to  the  fitted  dipole  model. 

The  principal  testing  carried  out  so  far  and  described  herein 
was  for  a set  of  reference  points  at  200-km  altitude,  distributed  over 
-60°(30‘’)60'^  in  north  latitude  and  0^(60°)300°  in  east  longitude.  For 
each  of  these,  a set  of  test  points  at  200-km  altitude  was  specified  with 
offsets  in  latitude  of  +5°,  0°,  -5°  and  in  longitude  of  +10°,  0°,  -10° 

(a  total  of  nine  test  points  for  each  reference  point).  A field- line  inter- 
section altitude  of  60-km  altitude  was  called  out,  in  separate  runs  for 
both  the  near  and  far  magnetic  hemisphere.  (One  additional  reference 
point  with  a set  of  test  points  near  the  south  magnetic  pole  was  also  run.) 
A check  on  the  field- line  intersection  locations  was  provided  by  separate 
runs  of  the  LINTRA  routine. 
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The  results  of  the  field- strength  comparisons  are  illustrated 
in  Fig.  7 in  the  form  of  a histogram  of  the  distribution  of  errors.  It 
will  be  seen  that  the  standard  deviation  is  of  the  order  of  1-2  percent. 
However,  some  test  locations  were  found  where  the  error  was  con- 
siderably larger  than  this.  The  geomagnetic  field  has  considerable 
deviations  from  a dipole  in  some  parts  of  the  world.  Nevertheless,  it 
is  felt  that  the  fitted  dipole  model  is  of  acceptable  accuracy. 


The  results  for  the  tests  of  field  line  intersection  locations  can 


be  summarized  by  stating  that,  for  intersection  locations  in  the  near 
magnetic  hemisphere,  the  average  latitude  error  was  0. 038°  and  that 
in  longitude,  0.019°.  However,  the  median  errors  in  both  latitude  and 
longitude  were  about  0. 01°,  showing  again  that  occasional  errors  much 
larger  than  the  average  occur. 


As  for  the  location  of  intersection  points  in  the  opposite  mag- 
netic hemisphere,  the  less  said  the  better,  in  general.  The  present 
ambient  geomagnetic  field  model  is  a local  best  fit,  and  that  is  not  a 
procedure  that  gives  a good  fit  in  the  large. 


wm 


Percent  Error 


Fig.  7.  Frequency  i)istributi on  of  Errors 
in  Total  Field  Strength. 
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APPENDIX 


In  this  appendix  are  included  certain  materials  of  interest 
only  to  those  who  wish  to  exercise  this  model  on  their  own  computer, 
and  who  presumably  have  a Fortran  card  deck  or  tape  available. 

Table  A1  contains  a definition  of  the  variables  used  in  the 
equations  of  the  text,  and  a cross  reference  with  the  Foilran  names  of 
the  variables  used  in  the  listing. 

Tables  A2  through  A6  contain  lists  of  the  input/output  quanti- 
ties for  subroutines  ONEMG5,  MAGFIT,  BFIELD,  CONJUG,  and 
MAGDRV. 

Tabic  A7  contains  a compile-and-run  listing  of  the  whole 
module,  together  with  the  input  and  output  for  the  test  problem  described 
in  the  text  of  the  report. 

Finally,  Table  A3  contains  a summary  of  our  experience  of 
the  running  times  of  the  various  routines  on  a CDC  7600  computer. 
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Symbol 

Fortran 

sin(\^) 

SINLTO  I 

cos(X^) 

COSLTO  ( 

‘^o 

PHIO 

M 

MUO 

X 

ANGS 

<P 

ANGE 

r 

RCUBE 

Br 

BR  \ 

Be 

BTHETA  > 

BPHI  * 

X 

CHI 

cos(x) 

COSCHI 

a 

ALPHA 

b2 

B2SQ 

I 

DIPANG 

D 

DECANG 

B 

BVAL 

Be 

RE 

^1 

AXATl 

Table  Al.  Symbols  and  Their  Fortran  Names. 
Drtran  Definition 


Sine  and  cosine  of  north  latitude  of  magnetic 
dipole  moment 


East  longitude  of  magnetic  dipole  moment 
Magnetic  dipole  moment 
North  latitude  of  field  point 
East  longitude  of  field  point 
Geocentric  radius  of  field  point 


Geocentric  spherical  field  vector  components 
(Bf,  positive  outward;  89,  positive  southward; 
and  B(p»  positive  eastward)  of 
IGRF  1965. 0 field. 


Angle  between  the  magnetic  dipole  moment 
vector  and  field  point 


Cosine  of  angle  between  the  magnetic  dipole 
moment  vector  and  field  point 


Angle  between  the  magnetic  dipole  moment  and 
geographic  north  pole 


Square  of  the  angular  component  of  the  magnetic 
field 


Magnetic  dip  angle  at  field  point 
Magnetic  declination  angle  at  field  point 
Magnetic  field  strength  at  field  point 
Radius  of  earth 
North  latitude  of  point  1 


Table  Al.  (Continued) 


Symbol 

Fortran 

Definition 

^1 

ALONl 

East  longitude  of  point  1 

^1 

ALTl 

Altitude  of  point  1 

^2 

ALT2 

Altitude  of  point  2 

^2 

ALAT2 

North  latitude  of  point  2 

^^2 

ALON2 

East  longitude  of  point  2 

AJUG 

Flag  controlling  which  magnetic  hemisphere  the 
location  of  the  intersection  point  is  calculated 

sin(Xj) 

SINZl 

Sine  of  angle  between  the  magnetic  dipole  moment 
vector  and  point  1 

COS(Xj) 

COSZl 

Cosine  of  Xj 

sin(X2) 

SINZ2  j 

Sine  and  cosine  of  angle  between  the  magnetic 

cosCXg) 

COSZ2  ) 

dipole  moment  vector  and  point  2 

sin(j^'j) 

SINPSI 

Sine  and  cosine  of  angle  between  the  field  line 
and  line  joining  the  magnetic  dipole  moment 

cos(^^) 

COSPSI 

with  the  north  geographic  pole 
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Table  A2.  0NEMG5  Subroutine  Input/Output. 


INPUT  VARIABLES 
TM 
RKM 
ST 
CT 
SPH 
CPH 

OUTPUT  VARIABLES 
BR 

BTHETA 

BPHI 

B 


Time  in  years  for  desired  field 
Geocentric  distance  of  point  (km) 

Sine  of  (geocentric)  colatitude  of  point 
Cosine  of  (geocentric)  colatitude  of  point 
Sine  of  (geocentric)  east  longitude  of  point 
Cosine  of  (geocentric)  east  longitude  of  point 


Radial  field  component  (gauss) 
Positive- south  field  component  (gauss) 
Positive-east  field  component  (gauss) 
Total  field  magmtude  (gauss) 
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Table  A3.  MAGFIT  Subroutine  Input/Output 


INPUT  VARIABLES 
Argument  List 


ALATF 

North  latitude  of  specified  point  P (radians) 

ALONF 

East  longitude  of  specified  point  P (radians) 

ALTF 

Altitude  of  specified  point  P (km) 

TM 

Time  for  desired  field  (years) 

OUTPUT  VARIABLES 

MAGLNK  Common 

MUO 

Magnetic  dipole  moment  (gauss  km^) 

COSLTO 

Cosine  of  north  latitude  of  magnetic  dipole 
moment 

SINLTO 

Sine  of  north  latitude  of  magnetic  dipole 
moment 

PHIO 

East  longitude  of  magnetic  dipole  moment 
(radians) 
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Table  A4.  BFIELD  Subroutine  Input/CXitput 


INPUT  VARIABLES 


Argument  List 
ANGS 


ANGE 


ALT 


North  latitude  of  field  point  (radians) 
East  longitude  of  field  point  (radians) 
Altitude  of  field  point  (km) 


MAGLNK  Common 


MUO 


COSLTO 


O 

Magnetic  dipole  moment  (gauss  km"*) 

Cosine  of  north  latitude  of  magnetic  dipole 
moment 


SINLTO  Sine  of  north  latitude  of  magnetic  dipole 

moment 


PHIO 


East  longitude  of  magnetic  dipole  moment 
(radians) 


OUTPUT  VARIABLES 


Argument  List 


BVAL 


DIPANG 


DECANG 


COSCHI 


Magnetic  dipole  field  strength  at  point  (gauss) 

Dip  angle  of  the  magnetic  dipole  field  at  point 
(radians) 

Declination  angle  of  the  magnetic  dipole  field 
at  point  (radiims) 

Cosine  of  angle  between  the  magnetic  dipole 
moment  vector  and  field  point 


g 


I 

A 

& 

I 

li 


I 

I 

m 

t 

I 


e 


Table  A5.  CONJUG  Subroutine  Input/Output 


INPUT  VARIABLES 

Argument  List 
ALATl 
ALONl 
ALTl 
ALT2 
AJUG 


MAGLNK  Common 
MUO 
COSLTO 

SINLTO 

PHIO 

OUTPUT  VARIABLES 

ALAT2 

ALON2 

S12 

BEBl 


North  latitude  of  point  1 (radians) 

East  longitude  of  point  1 (radians) 

Altitude  of  point  1 (km) 

Altitude  of  point  2 (km) 

1.  *•  Calculates  latitude  and  longitude  of 
point  2 in  same  magnetic  hemisphere 

-1.  - Calculates  latitude  and  longitude  of 

point  2 in  opposite  magnetic  hemisphere 

3 

Magnetic  dipole  moment  (gauss  km  ) 

Cosine  of  north  latitude  of  magnetic  dipole 
moment 

Sine  of  north  latitude  of  magnetic  dipole 
moment 

East  longitude  of  magnetic  dipole  moment 
(radians) 


North  latitude  of  point  2 (radians) 

East  longitude  of  point  2 (radians) 

Path  length  along  the  field  line  from  point  1 
to  point  2 (in  units  of  the  equatorial  radius  of 
the  traced  field  line) 

Ratio  of  the  equatorial  value  of  the  field  to  that 
at  point  1 for  the  traced  field  line 
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Table  A6.  MAGDRV  Input  Quantities  — START  Namelist 


ALATFI  Array  of  north  latitudes  of  fit  points  (deg) 

ALONFI  Array  of  east  longitudes  of  fit  points  (deg) 

ALTFI  Array  of  altitudes  of  fit  points  (km) 

NFIT  Number  of  fit  points 

TM  Time  at  which  to  evaluate  exact  field  (years) 

RLATS*  Array  of  north- latitude  deltas  of  test  points  (deg) 

RLONS*  Array  of  eust-longitude  deltas  of  test  points  (deg) 

RALTS*  Array  of  altitude  deltas  of  test  points  (km) 

NRS  Number  of  test  points  relative  to  a fit  point 

RCONS  Array  of  test  altitudes  for  intersection  calculations  (km) 

AJUGS  Array  of  calculation  options  for  conjugate- region  inter- 
section calculations: 

1.  - Calculate  intersection  point  in  same  magnetic 
hemisphere. 

-1.  - Calculate  intersection  point  in  opposite  magnetic 
hemisphere. 

lOPT  MAGDRV  calculation  options: 

1 • Calculate  only  magnetic  dipole  field  at  test  points. 

2 - Also  calculate  location  of  intersection  points. 

3 - Also  calculate  magnetic  multipole  field  at  test  points. 

*The  (input)  locations  of  the  test  points  are  relative  to  the  fit  point. 
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Table  A7. 


Compile-and-Run  Listing  of  the  Ambient  Magnetic-Field 
Module,  with  Input  and  Output  of  Test  Problems. 


MAGDRV 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


M*(;n«V(TNPliT,nuTPUT,T*PES*H»’UT,  UPFbsil'ITPUT) 


THIS  PBOGRAM  FkEhCISFS  TmE  4mBIE»,T  haShETIC  EIELU  KOOEt,  THE 
HflDFL  CONSISTS  OF  fOijB  PtlOTtNES,  h.agPIT  FITS  i DlPOtf  FIfLO  TO 
the  HAGNETir  FIFLO  AT  A GIVEN  POINT, WHICH  SMIJULO  BE  NfA» 
the  center  of  1 -IE  HEf.I'JN  OF  INTEREST,  THE  PEVACT*  FIELD  AT  THE 
POINT  IS  rAlt'iLATFO  FROH  ONfHGS#  A MJt>EL  OF  The  INTERNATIONAL 

gfohagnfttc  reference  fiflo.  epoch  IRES.O, 

(SEE  fll  E.G,  STaSSINOPOIjLlIS  AND  R,D,  HEAD#  NASA  REPORT 
NSSnc  7?-12,  ALI.ma6»«DALhG.LINTRa  CORPOTEH  PROGRAHS  EOR 
GFiihaGmFTTC  FtElO  AND  EIELD-LTnF  CALCOLATtONS#  EEBRijARy  IR72 
AND  (?)  J,c.  CaTn  and  S,J*  CaTn,  NASA  TN  D»62)7,  derivation 
OF  THE  TnteRnatTONaL  GEOhaGNETIC  BLEERENCE  field  (lGREtlO/6E)>» 
All(.UST  1971,) 

A third  RouTiNj,  reield.  calculates  the  magnetic  field  strength 
FOR  any  given  point  for  the  fitted  DlPOLt,  ThE  FUURTm  ROUTINE 
calculates  The  location  of  a point  with  a given  altitude  hh1C’< 

IS  ON  THE  SAME  FIELD  LINE  AS  SOHE  specified  POINT  FOR  THE  FITTED 
OTPOLF  field. 


Input  p 

ALATE 

Al  one 

ALTFI 
NF  it 

fH 

»LATS 

RI.UNS 

halts 

NRS 
PC  ON  3 
A.iur.s 


T . 


ARAhETERS  (NAhFlIST  START) 

I  . ARPaV  of  north  LAtlTUDES  OE  FIT  POINTS  (DEG) 

array  of  east  longitudes  of  FIT  POINTS  (DEG) 

array  of  ALTITUI'CS  of  fit  points  C^M) 

NuHrFH  of  fit  points 

TlHg  aT  hhICH  VO  evaluate  exact  FIELD  (TEARS) 
ARRaV  of  north  LATITUoFS  CjF  test  points  (DEG) 

array  of  east  LUNgITIIdfS  of  test  points  (DFG) 

ARRaT  of  AlTITOues  OF  TEST  POINTS  (RH) 

NOHrtFR  OF  TEST  POINTS 

AHRaV  of  TFST  field  line  altitudes  (KM> 

ARPay  of  test  field  line  calculation  options 
t,  calculates  InTeRSFCTIOn  point  in  same 
Magnetic  hemtspepe 

-1,  calculates  intersection  in  ('PPuSITE 
maGNETTC  hEMISPePF 


CAUTION  - location  of  opposite  hemisphere 
TnTeRSFcTIUnS  may  not  he  accurate 


inPT 


calculation  Options 

1 - calculate  ONLV  dipole  » field  aT  TEST  POINTS 

2 • also  calculate  location  of  Intersection  Ptu^ts 

3 - also  calculate  HulTIPOLF  b field  at  test  points 


FIT  dipole  to  point 


REAL  M||0 

common  /maglnk/  MUO»r.OSLTO»8INLTo,PHlO 
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NEWM4G.2 

NEw^AG, 5 

NEmMaG.R 

NEnHAG.S 

nEhmaG,^ 

nEmmag.t 

NEmHAC.B 

NEhMaG.R 

NEhHaG.10 

NEmMaC,11 

nEhMaG,12 

NEwMAG.li 

NEwMAGtIU 

nEnMAG,)5 

NEwHAG,16 

NEi*maG,17 

NEWMAG.ia 

nEnmaG.IR 

NEhMaG.EO 

NERMaC,?! 

NFwMAG,?? 

NEmMA6,25 

n£nRAG,2U 

NEhhaG,25 

NEMHA6,2b 

nEmHaC,27 

NEMRAC,2a 

NE»(HaG,29 

nEnRA6,50 

NE"MA6,SI 

NErMaC,32 

NEhhaGjSJ 

NEnMAG,3<I 

nEnHaG.JS 

NERMAG.Sb 

nE».haC,37 

NEhMAG.jS 

N£*.HaG,J9 

NERMAG.ao 

NErMaC.RI 

nFnMa6,U2 

NEnma6,«3 

NErMAG,<IU 

NEwMaG,U5 

NENMAG,a6 

NEhMa6,<I7 

NEHMAC.Oe 

NE!i<MaG,«9 

NEmMAGoSO 

nErmag.SI 

nE<iMaG,52 

nErRAG.SJ 

NEpMAGfSR 

NE"MAG,5S 

NERMAG,Sb 

Mhma6,S7 


MAGDRV  (Cont'd) 


fi 

n 


11 

20 

21 

^^ 

u1 


hi 


COMMllN/CNStNT/  «tCM»1'T»M*LFPI»Fnil»PI,r.rt*V2(?) 


NEwM*6,S8 
ht»»M*G,S9 
NEwMiC.fcO 
NEMMAC,6i 
NEi»HAG,62 
NEMMAG,6i 

NAMFLfST  /ST*aT/*L*TFT,*LONFl»*LTFi.NFlT,T-<,HL*TS.«LON8,»*L18»NRS,HEi»M*S,6a 

t.C 


OTMENSiriN 

OTMERSION 

dimension 


Rl  ATsrsO)  .RLONSCSfl)  >iiCON8(SO)  ,*JUGS(S0) 
»*LTS(S0) 

»liTFT(SO)i*LCINFI(SO)»*LrElfSO) 


»RCONS,A.IUGS.  ICIPT 


NE  * 
RADS 


1 .nE»0S*»fcCM 
X MALfPl/90, 


a? 

«•) 


READ  IN  Data 


BFAD(S»STANY) 
wRITE(6. START) 


LOOP  OVER  fit  PijTnTS 


on  POO  jj»i,nfit 

Al  ATF  s At  ATFI(JJ) 

ALONF  3 ALONfT(JJ) 

ALTF  s AUTFt(Jj) 

wRITEIo. t 000) ala TF, ALONF, AUTF»TM 
Al  ATF  » ALATF'RAOS 
ALONF  * AlONF*RAOS 

000  FnRMAT(1Hl,1}H  location  UE  PUInY 
T.FR.J.ftH  (0FG),/,13H  longitude  * 

»»Ffl,?,6H  (0FG),/#13H  time  • »F«,2,6M 


THAI  IS  FITTED»//»12H  LATITUDE 
,E8,2,8M  (0EG),/,13H  ALTITUDE 
LVR8)»///) 


FIT  dipole  to  riven  field  point 


* ,E13tt>> 


1 


hS 

70 


70 

77 


103 


call  MAGFlTf ALaTF, ALONF, ALTF,Tm) 
*iRlTtfb.lOOl)M||0,COSLTO.SlNLTO,PHIO 

001  FORMAT!//, ^SH  fitted  niPUte  parameters,//# lOH  MUO 
»13H  GAUSS  '<M**S,/,10H  Cn8LT0  * ,€  1 3,h, /,  1 OH  SlNLTO  ■ ,Et3.6,/, 

»ioM  PHio  s ,E 3,t),a(iH  longitude  East  (Raoians),///) 

10  oRITE(6,1002) 

002  FnHMAT{lH*,lX,'»H  TEST  L4T,2X,9h  TEST  L0N,1X,9h  TEST  ALT,2X, 

$ Oh  OTPOlE  8,1x,7H  0IPANC,5k,7h  OECAnG) 

IFIIUPT  ,l  T.  2)  GO  TO  50 
?0  iiiRITE(6,1003) 

003  FORMAT!  IMF, 62X,  1 OHINtFRS  *LT,2x,6H  AJuG  ,«X,  IOMINTERS  LAT»21<» 
*10mINTER8  LON) 

IF  IinPT  ,LT,  3)G0  to  50 
WR1TE!6,100«) 

FnHMAT!10'ix,7W  IGRF  H,2X,SH  PERCENT) 

continue 


30 

ona 

so 


LOnP  OVER  test  points 


103 


DO  <100  J3l  ,NRS 


fALC'iLAIF  DIPOLE  FT£LO  VALUE  *T  TEST  POINT 


lOh 


ANGS  a AlATF  ♦ RlATS!  I)»hAOS 
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nEnm«c,65 
NEhHAG|66 
NEmMaG,67 
NEnMaG,68 
NEpHaG,69 
NEnMaG,70 
NEnHAG.TI 
NEwHaG.72 
NEmMaG.73 
NEhMaG.74 
NEpmaG,75 
NFRMAG,76 
NE**maG,77 
NEmmaG,78 
NEmMaG,79 
nEnMAG,S0 
NEnMaG.SI 
nEhHAG.S2 
NEmhaG,S3 
sNEwma6,BU 
« NEnMaG,65 
NENNAR,6b 
NEmMaG,87 
NEnNAG,S8 
nE»>maG,89 
nEnMaG,90 
NEnMaG.RI 
nEnMaG,92 
NENMAG,9i 
NE*<MaG,9U 
NEnMAG,95 
NEmMaG,96 
NEmMaG,97 
NEwmaG,98 
NEmMaG,99 
NEhMaG, too 
NEwMaG.IOI 
NEwMaG,102 
NENHA6,tOS 
NFwmaG.IOO 
NEnMaG,105 
NEmmaG, lOb 
NEMMAG,t07 
NEmMaG.IOB 
nEwMaG,109 
NEmmaG, 1 10 
NEmHaG.UI 
NEnMaG.I 12 
NEwMAG.in 


MAGDRV  (Cont’d) 


tin  ANGF  s alHNF  ♦ PunNiiC  ))*(;*0S 

113  *L1  * *LTF  ♦ R*LT9(J) 

llS  CALL  BF!£LO(ANGS|ANGe«*lT*B«OTPANC*nLCANG.COSCH!) 

120  ANGSO  « ANGS  / RAOS 

t2A  AMGFO  a ANCF  / RADS 

127  DtPAhf,  « OIPaNg  / PADS 

131  OFCANG  a PECANg  / HAOS 

132  MRlTE(b«1005)  AMGSO«ANGEO*ALT*H.DIPAnG«PECAng 

1005  F0HMAT(1M*,1(?);,FP,2),2X,Fe,5»3)f»F7,2,2if.Fe,2) 

152  IFdOPT  ,LE.  nr.O  TO  OOO 

C 

C CAlCULATF  LOCATirSN  OP  POINTS  at  A^TITUPE  RCONS  nHJCM  are  on 
C FIELP  LINE  AS  the  test  POINT 
C 

15A  Ajur.  a AJUCS(J) 

157  ALT2aRC0NS(J) 

161  CALL  CaNJU6(ANGS,ANCE,ALT.ALT2«AJiiG.ALA2«ALD2,si2,BFB1) 

171  ALA2aALA2/RADS 

173  AL0?aAL02/RARS 

170  NRITE(6«ia06)ALT2»AJUG»ALA2,ALO2 

1006  FOHHAT(lHA«60X,F10,2«F9.2*tt)(«2(Fin,3>2k)) 

207  IFdtlPT  ,LE.  2)  CO  TO  oOO 

C 

C calculate  aEXACT*  FIElO  PROh  ir.RF(196bt)  AT  TEST  POINT 
C 

212  RKM  a RF  ♦ alt 

213  COLAT  a HALFPI.ANGS 

216  ST  a SIN(COLAT) 

217  CT  a COSICOLAT) 

221  SPH  a SINIANGF) 

223  CPH  a COS(aNGE) 

226  CALL  0NEMC5dM,RaM,8T»Ct.SPM,CPH,BK,BT,LP,BExACT) 

23T  pelt  a(R  • REXACTI/BEXACTaIOO, 

203  mRITE(6.1007)BEXACT»OELT 

1007  FORNAT(103x,2X,Fi0,3,1X»F10,3) 

252  000  continue 

253  900  CONTINUE 

257  t<RlTE(6,100B) 

lOOB  FORm*t(///,20H  End  of  test  PROBLEM) 

263  ST(jP 

265  END 


BFIELD 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


1? 

l« 


SDPRUUTINE  RFTELO(ANGn,*NGt»ALT«BVAL>DIPANa,Di;C*MG,CUSCHn 


TMIS  routine  calculates  THfc  AMBIENT  MAGNETIC  FIElO  at  a POINT 
FROM  THE  m*GnTt"OF  AnO  OlRtCTlOM  OE  THE  MAGNETIC  DIPOLfc  MOMENT, 
AND  THE  location  of  the.  POINT 


INPUT  parameters 
ARGUMENT  LIST  • 

*ng-  . north  Latitude  of  field  point  (radianbj 
ANGF  East  longitude  of  field  point  (RAOlALST 
ALT  - altitude  of  field  POINT  (KM) 


NEmmaC.IES 

NEwmaG.ISE 

nEwmaG.IST 

NEhMaG.ISP 

nEmmag.ISR 

NEmMaC, tbO 
NEMMAG.lbt 
NEHMAG.IfrZ 
NEmHAG,163 
NEMHAG.tbU 
NEmmaG, J65 
NENMAG.tbb 
nEnHaG, 16T 
NF.MMAG.lbS 

,170 


maglnk  common 

MUO  « magnetic  OIpOlF  moment  (GAUSS«KmT) 
cosLTo  - rnsiNE  of  north  latitude  of  magnetic  dipole  moment 

SINLTft  • Sine  of  north  latitude  of  magnetic  dipole  moment 

pHio  • East  longitude  of  magnetic  dtpole  moment  (radians) 

CNSTNT  common 

PI  • S.luibRabSlSB’R 

RFCM  • Earth  Radius  (cm) 


nEhmaG 

NEhmaG 

nEmmaG 


NEmmag,171 
NEmMaG,172 
NE»>maG,175 


NEhmaG,  t7<l 
NEnMAG.t  7S 


NEhMaG.Wd 
NEhmaG,  1 )'7 

Aif  A ^ A « A 


OUTPUT  parameters 
ARGUMENT  LIST 


BVAL 

OIPANG 


Magnetic  field  strengmt  at  field  point  (gauss) 

DIP  ANGLE  OF  magnetic  FIELD  AT  FIELD  POINT  (RAD), 

DIP  (OR  INCLINATION)  IS  THE  VERTICAL  ANGLE  MEASURED 
FROM  THE  mORWONTaL  aT  ANY  POINT  TO  THE  (VECTOR) 

LINE  OF  force  Through  that  point,  it  is  positive  in 
The  northern  magnetic  hemisphere  and  negative  in 
The  southern  magnetic  hemisphere. 

Declination  of  magnetic  field  at  field  point  (rad), 
declination  (OR  Variation),  the  angle  bethfen  the 
geographic  and  magnetic  meridians  at  a point,  is 
Positive  if  the  compass  nefdle  points  east  of 
geographic  north, 

COSCMI  • cosine  of  MaGneT1c*DTPOLF  COLaTITUOE, 

COSCHI  IS  negative  in  The  northern  magnetic 
hemisphere  aNq  positive  in  the  southern  magnetic 
hemisphere. 


CFCANG 


real  MUO 


COMMON  /maglnk/  mU9»COSLTO,SINlTO,PMIO 

CnMM0N/CNSTNT/RECM,PI,HALFPT,FoUHPl,6RAVZ,GZHE2,B0LTZK,GAMj,GMll 

,PMNIT,PMOxV 


HF«RECH*1 ,£-S 

calculate  SINE  and  COSINE  OF  NORTH  LATITUDE  OF  HELD  POINT 


COSLAP  B COS(AnGS) 
STNLAP  a SIN(AnGS) 


• 1 r f 

NE«MaG,178 
NE*>maG,I7R 
NEhHaG  . < .->C 
nEmmaD,  •.01 
nEmma.-;,  ii-? 
nEm'^iaG, 
NEHf  ^i;,  ISO 

NE> -«5,ies 
nEwhaG,166 
NEhmaG, 167 
NENMAG,lBe 
NEhmaG, 169 
NEhmaG, 190 
NEhmaG, 191 
NEhmaG, 192 
NEhmaG, 195 
nEhMaG,190 
NEHMAG, 19G 
NEHHAG,t96 
nEhmaG,197 
NEhmaG, 19S 
NEhmaG, 199 
NEhmaG, 200 
NEhMaG,201 
NEhmaG, 202 
NEhmaG, 205 
NEhmaG, 200 
NEnMA6,20S 
NEhmaG, 20(1 
NEhmaG, 207 
nEhmaG,20S 
NEhmaG, 209 
NEhmaG, 210 
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BFIELD  (Cont'd) 


?t  OEULUN  s ANSE  • PHIO 

C 

C C*LCI.'UATf  8lM€  AKO  COSlinf  Of  AM6LE  BETWEEN  HASNETIC 
C 0I?0LE  mOhEMT  AmO  flELO  BOInT 

27  CnSCMI  « COSLAp*cnSl.TO*Co8CAhCE  • fHloj  ♦ SXULAf aSINuTO 

at  5tl»CWI  « SaOTtt.  • CUSCMIaCOSChI) 

c 

C rALCliLATE  CUBE  Qf  GEOCCnTRXC  RAOluS  Of  POINT 
C 

St  RCUBE  ■ «E  ♦ AI.T 

S3  RCU8E  « ReuNE*RCuSE*Reu8E 

c 

C total  flEL"  strength 

c 

Sa  BTEBM  « SQRT(S.*COSCHr*CoSCMI  ♦ l,) 

(>S  BVAL  » nuO/RCUSEaBTERN 

c 

C haSnETTC  OIR  ANSLE  at  POINT 

c 

*7  StMOIP  * ?,  * COSCMI  / BiERN 

72  BIPANC  * .ASINtSlNBlPJ 

C 

C magnetic  OECUInaTION  ancle  at  point 

c 

77  STnPSI  p COSlTo  • SlNtOEtLON)  f SINCHI 

t03  OECANC  « ASlNlsXNPSt) 

1t5  COSPSI  « -SINLTO  ♦ COSCNI  • SINLAP 

120  If CCOSPSI.lt, 0.)  OECANC  * SI6N(PI,SINP$I)  • OECANC 

C 

125  RETURN 

t2A  END 


NEpMac.Zlt 

NEwM*c.212 

NEwNAC.213 

NEwNAC.ZlA 

NEwNAC.2t5 

NEANAC.2t6 

nEmMaC.ZIT 

NE«MAC,2te 

NEWPAC.21< 

NEMMaC.220 

nEwMaC.221 

NEnnaC,222 

NEwMAe.223 

NEwMaC,22A 

NEwnaC.225 

NEMMAC.22B 

».rwHA6,227 

NEwmaC.228 

NEpnaC.22* 

NE»NAe.230 

NEpnaC.251 

NEhnaC.232 

NEWMAC.23S 

NEHMAC.25A 

nEppaC.235 

NEPNAC.23B 

NEWNAG.23T 

NEWNAC.23S 

NEIINA6.23* 

NEWMAC.2A0 

NENMAC.2«t 

NEhnaC,2a2 


% 
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R 

'M 


BLOCKH 


BLOCK  D«T*  BLOCKH 


ImXT1*uI2C  named  common  CONSTANTS  AND  DEFAULT  VALUES 


SET  0^  CHEm  quantities  ApPEARIn£  UNDER  VARIOUS  CHEM  QRTIONS 
DEPOSITION  COEEEICIENTS  fOR  PROMPq 
COMMON/DEPOAT/sTCU(5,*),R»»Eff t5»«)»£R6ut5),TMRE8M(a),8PlNTM, 

1 SPIRT, *MU1 

mathematical  and  sEOPMrsiCAL  constants 

C0MMUN/CNSTNT/RE,Rl*HALEPI»r0URRl,CRAVZ,&ERE2*S0LTZK,SAMt,6MlI 
1 ,Rmn1T,omOxT 


CNSTNT 


DATA  RE  /6, 

1 PI  /3. 

3 halfPi/u' 

5 EOljRPl/12 

7 CRAvZ/RBO 

5 ROLTZK/i, 

? PMNJT  /2, 

3 FmOxV  /^,i 

h GaM|/o.5/ 

OEPDAT 

DATA  SI6U  /O.O 

1 2.0E«tB 

2 0.0 

3 0,0 

« ER6u  /t.TftZE* 
$ TmRE8m/|,s*«E 
fc  SPlNTM/i.OOEl 
END 


3*7650E^0«/, 

I*1SR2653S8P5/, 

S7079632s79aR/, 

•566370*103992/, 

,665/* 

3B05«E*16/, 

929763E-2S/# 

656RSeE*23/* 


,B.0E*l8,2,0E*t7,2,5E»i7,2,0E«17, 

,l,6Eol7,2.0E»l7,2.5E.l7,2.0E»l7, 

,t.0E*l7,l,0E»l7,t.0E«l7,t .0E»17, 

»3.2Emi8,5,2E.1B,9,OE*18,9,oE«IB/, 

lt,2*37lE«|l,2,S95E.il,3*525E*il,5,767E*ll/, 

•lt,B,202E*12,2.33lE»n,2,lB2E-U/, 

l/,8PlNT/2,07Elt/ 


nEmmaC.ZOS 

NEMHA6,2a6 

NEmmaG,205 

NEHMA6,2a6 

NEnHa6,247 

nEmMa6,2«8 

NE«MA6,2a9 

NEMMA6.250 

NEmma6,2S1 

NEmma6,2S2 

NEMMA6.253 

NEmMa6,?5<i 

NEMMA6.255 

NEmMaC,256 

nEmM46,257 

NEmmaC.ZSB 

NEmma6,?S9 

NEmMaG.260 

NEMMAC,26t 

NEMMA6.262 

NEMMAS.26S 

NEnmaG.26R 

NEmma6,265 

NEmmaC,266 

NE*ma6,267 

NEmma6,26S 

NC«ma6,269 

NEmMa6,270 

NEmmaG,271 

NE»mac,272 

NEmmaC,273 


CONJUG 

SUBROUTTNF  CONjuG(A'kAT|»Al.UNi.«Lt  t **LT2«*Ju6,*tAT?>  AL0N?,Sl2fREBnNEMHAG,3ra 

C NEnmaC,27S 

C********  nCmMAG,276 

C NENHiG,277 

C THE  BOUTIHE  Cf’NJOG  CALCULATES#  Eo#  A GIWEoi  ALTITUDE#  THE  lOCATIOH  NEAMAG,27e 

C (LAT,LON)  OF  THE  BOlNT(2)#  WHICH  TS  ON  THE  SAME  MAGNETIC  DIPOLE  NE><maG,?79 

C FTELO  LINE  AS  S^^HE  0THE»  GIVEN  POTNT(!),  IT  ALSO  CALCULATES  THE  NEMMAG,2eo 

C FIELO^LINE  DISTANCE  eFT-EEN  POlNT(I)  AND  POIn7(2)#  IN  UNITS  OF  NEii<MAG,2ei 

c the  equatorial  Radius  of  the  traced  line#  and  the  ratio  of  the  nemmag,282 

c equatorial  Value  of  the  field  to  that  at  puintiu  for  the  traced  NENMA6,2e3 

C line,  nEnmaS,2«R 

c NEMMAG.aas 

Ca******a  NEmhaC,286 

C nEhMaG,287 

C INPUT  parameters  nEhMaG,288 

c argument  list  • NEnMAG.289 

c ALATt  - NORTH  latitude  OF  POINT  • (RADIANS)  nE»«maC,290 

C ALONl  • East  longitude  QF  point  I (radians)  NEmMaG,291 

c ALT1  • altitude  of  point  1 (AM)  nEmmaG,292 

C ALT2  • altitude  OF  POINT  2 (KM)  NEmMaG,293 

C AJUG  • NEhMaG.29« 

C 1,  calculates  LOCATTUN  (LAT.LON)  of  point  2 nEwmaG,295 

C IN  SAME  M46NfTlc  hemisphere  NEi»MAG,296 

C -1,  calculates  location  (LAT,L0N)  of  point  2 NEnMaG,297 

c IN  OPPOSITE  Magnetic  hemisphere  NE»iMAG,298 

c nE»MaG,299 

C A«**«*******«*«**A*******A*««««**»*A*************  NEmHAG,300 

C * • nEomaG.SOI 

C * CAUTION  • LOCATION  UF  OPPOSITE  HEhIsPHERF  * nEnmaG,302 

c A Intersections  mat  not  be  accurate  * nEmmag.sos 

c * * NEHMAG.30R 

c A«************A***A*******A**********************  nEnmaG,30S 

c nEhhaG,308 

c MAGLNK  common  nEixmaGiSOT 

C HUO  • magnetic  OIpOlE  moment  (GAU5S*km3)  NEmMAG,30B 

C COSLTO  ■#  COSINE  OF  NORTH  LATITUDE  OF  MAGNETIC  DIPOLE  MOMENT  nE#iMaG,309 

C SINLTO  V SINE  OF  north  LATITUDE  OF  MAGNETIC  DIPOLE  MOMENT  NEi»maG,310 

C PHIO  - East  longitude  of  magnetic  dipole  moment  (RaOIAnS)  nEmhaG.SH 

C CNSTnT  common  nE"Ma6,3I2 

c RECM  • Earth  Radius  (cm)  ne«mag,5i3 

C PI  • 3, iai 5926535898  nE»maG,310 

C nEhmaG,315 

C OUTPUT  parameters  nE«<maG,316 

C ALaT^  - NORTH  latitude  of  POINT  2 (RADIANS)  NEmhaC.SIT 

C AL0N2  • east  longitude  OF  POINT  2 (RADIANS)  NE«MAG,3t8 

C S12  • distance  along  field  between  point  I AND  POINT  2 (IN  NEwHaC,S19 

C UNITS  OF  the  equatorial  value  of  The  traced  LINE)  NEnHAG,320 

C PFei  • RATIO  OF  THE  equatorial  VALUE  OF  THE  FIELD  TO  THAT  AT  nEbMaG,321 

C POINT  1 FOR  THE  traced  LINE  nEwmaG,322 

C NEHHA6.323 

COH-ON  /MAGLNK/  mU0.COSLT0#SINlT0.PmI0  NE"Ma6,32« 

COMMON/CNSTNT/  RECM,PT,HALFPl#FnuRPI#GRAV2(r)  NEwMaG,325 

C NEwhaG,326 

REAL  MOO  NE"MaG,327 

C NEmhaG,328 

RE  « l.nE«05*RFCH  NEwmaG,329 

la  TWOMI  s 2,'*PI  nEwmaG.SIO 

C NEWMAG,3il 

C calculate  sine  and  COSInE  of  Angle  between  magnetic  dipole  moment  nEwmaG,332 
C AND  POINT  1 nEhmaG,333 

C nEwmaC,33R 
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CONJUG  (Cont’d) 


f- 


IT 

20 

2S 

51 

41 

45 


4S 

55 

5T 


65 


72 

102 


i06 

11? 


117 

122 

12« 

127 

141 

145 

147 


c 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


C 

C 

C 


1 

153 

1. 

155 

15b 

s 

IbO 

1b4 

1 

171 

200 

213 

&- 

1 

kt. 

1 

217 

223 

223 

& 

A- 

COSl T1  • C0S{5l*T1) 

JTNLTI  • 8IN(*tiTl) 

DTFF  ■ ALONl  • I*HI0 

COSZl  • S1»<LT0*8INUT1  ♦ CO«L70*CO8t.Tl*COS(01FF) 

STNZ12  * i«  - coszi«ca$zi 
SfNZt  s SQfiT(S!NZ12) 

CALCUlATC  SIHE  AND  COSlAlC  OF  ANCLE  BET"EEN  HACnETXC  OlPQLt  HOHENT 
AND  POINT  2 

8INZ22  ■ 81NZ12*(»E  ♦ ALT2)/tPE  ♦ *UT1) 

8INZ2  ■ 8QP7(8INZ22) 

Cn8Z2  ■ 80«7(1,  - 8INZ22) 

OETERMTnF  sign  (IF  C082Z  • tF  POINT  2 18  IN  THE  SAhE 

HAGnFTIC  MEhISHEPE*  TMEn  SICN  OF  COSZ2  IS  THE  SAME  *8  THAT  OF 
COSZl  fCOSl'^E  OF  ANCtE  0ETHECH  MAGNETIC  OlPOtE  MOMENT  AND  POINT  1), 
IF  IN  OPPOSITE  magnetic  hEMiSPmfPe#  THEN  SIGN  OF  COSZ2  IS  OPPOSITE 
SIGN  OF  COSZl 

CnSZ2  a SI6N(CoSZ2«Aj(lG*COSZt) 

calculate  sine  *no  cosIne  of  angle  betmeen  field  lIne  and  line 
JOInTNG  magnetic  OIPOlE  HOHENT  MlTH  geographic  north  pole 

SINPSI  a COSlTI  a sIn(OIFF)  / SlNZl 

CnsPSi  a (sINLtI  • CO8ZlA8INLT01/tSlNZi*CO8LT0) 

CALCULATE  NORTH  LATITUDE  OF  POINT  2 

SINLT2  a COSZ2  * SINtTO  ♦ SInZ2  a COSLTO  * COSPSI 
ALAT2  8 ASINISlNtTZ) 

calculate  EAST  longitude  of  point  2 

SiNOIF  a SINZ2  * 8INFSI  / COSlALATZ) 

OTFF  a ASlN(SlNOlF) 

CnsSGN  a COSZ2  • 81NLT0*81NLT2 

IFICOSSGN  ,LT«  O.IOIFF  a SlGNtPf ,SlNOIF)  • OIFF 
AL0N2  a PMlO  ♦ DIFF 

IFULONZ  ,LT,  0.)  ALON2  a AL0N2  ♦ THOPI 
AL0N2  a AMOO(ALON?,TMnPIi 

calculate  dipole-field  path  length  BEThEEn  POINT  1 AND  POINT  2 (IN 
UNITS  OF  the  equatorial  RADIOS  OF  THE  TRACED  FIEuO  LINE) 


R5  a SQRTtl.) 

RSHui  a R5  * COSZl 
R3MU2  a R5  • COSZ? 

RlPRl?  a SORTll,  ♦ R3NU1  • 
R1PR22  a SQRTd.  ♦ R5mu2  • 


SlE  * ABS(R3M'Ji  A RtPRiZ  ♦ 
S2E  a ASSIRIMU?  • R1RR22  ♦ 
812  a ABSTSIE  • AJUC  * SzE) 


R5NOIJ 
R5HU2) 
AtOCtRsHUi 
ALOGlRTHUi! 
A R5/6. 


PlPRlZ)) 

RiPR22)> 


C 

C 

C 

C 


calculate  RATIO  OF  EQUATORIAL  VALUE  oF  THE  FIELD  TO  THAT  AT  POINT  1 

for  thf  tracfo  field  line 

BEBI  a SInzi2a*5/RjPrj2 

RETURN 

END 
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nEmm*G«555 

nEmmac.556 

NEmMfG,55T 

NEnmac,556 

nEmmaC,55R 

NEnHA6,540 

NEaHACt541 

nEmnaC,542 

nEmhaG,545 

NEMMAC,54a 

NENMAGt345 

NEmmaC.IUN 

NEmMaC.IuT 

NE«haG,348 

NEwmaG,549 

NEmMaC,350 

NEnHaC,351 

NEmmaC,352 

NEmmaG,353 

NEnma6,354 

NEnmaC.355 

nEmMaC.SSA 

NENMAC.35T 

nEmHaG.358 

NEnMa6.359 

nEaHaG,360 

nEwMaG,3B1 

nEmma6,3B2 

NENMAG,3k3 

NENMACt36a 

NEMMAC,3fc5 

nEmmaG»36B 

NEMMAGtSBT 

NCMMAG,3aB 

NEHMAG(3b9 

NEmmaG,370 

NEnmaGj371 

NE*tMAG,372 

NEwMaG,)73 

NEMMA6(374 

NEMMAG't3?5 

NEHMA6,S7b 

nEmma6,377 

NEMMAC.37B 

NEmma6,379 

NEmmaC.SBO 

NCmMa6,3B1 

NEMMAG.3B2 

NEmmaG,363 

NEMMAG.3B4 

NEMMA6.3B5 

NEMMA6,3Bb 

NEaMAG,3e7 

NEHMA6,3eB 

KEmha6,3B9 

NENMA6.390 

NE«|NAC,391 

nEmhaG.392 

nEhhaG«393 

nEmha6,594 

nEmhaG*395 

NC*»HA6,39b 


¥ 


MAGFIT 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 


T 

It 


1^ 

i« 

16 

20 

22 

2T 


c 

c 

c 


c 

c 

c 

c 


SUBROUTINE  HAOFIT ( aLaTE, alone, *LTF,TM) 

* • * * 

THIS  routinf  Fits  a otpole  field  to  the  local  magnetic 
field  at  a specified  POInT  p,  p i^  given  BV  ALATF,AlONF,ALTF, 

THE  magnetic  field  At  p is  found  from  model  S of  STASSINOPOULOS 
models,  model  S is  IGRF  10/6S,  reference  • STASSINQPOULOSf  E.C. 
AND  G.D.  head#  ALLMAC,  field-line  CALCULATION.  NaSA-GODDARD  SPACE 
FLIGHT  center.  NSSOC  T2-12.  FEBRUARV  1972, 


INPUT  parameters 

ALATF  • GEOCFNTRIC  NORTH  LATITUDE 

OF  sFECIFIED  POINT  P (RADIANS) 

ALONE  - GcOCENTRIC  EAST  LONGITUDE 

OF  SFECIFIEO  point  P (RADIANS) 

ALTF  • altitude  of  specified  point  P (KM) 

TH  • time  for  desired  field  (TEARS) 

CNSTNT  common 
HAlFPI  • PI/2 

PI  • 1.t«lS926SSS09a 

RFCM  • earth  Radius  (CH) 

OUTPUT  parameters  (To  maglnk  Common) 

MUO  • magnetic  dipole  moment  (CAUSS-KMS) 

COSLTO  • cosine  of  NQRtM  LATITUDE  OF  MAGNETIC  DIPOLE  MOMENT 

SINLTO  • SlNg  OF  north  latitude  OF  MAGNETIC  DIPOLE  MOMENT 

PHlO  • EAST  LONGITUDE  OF  MaGNETIC  DIPOLE  MOMENT  (RADIANS) 

rfal  lamoa.lamarg 

REAL  Mgo 

COMMON  /maglnk/  MU0»CO3LT0.SINlT0.PM10 

COMHON/CNSTNT/RFCM,PI.MAlFP1.FoURPI,CRAV2.GZRE2»BOLTZK,GAMI,GM1I 

1 .pmnit.pmdxt 

REpRECM*I,E-5 

TmOPI*2,*PI 

PIOVZaHALFPI 

SET  UP  INPUT  FOR  C*LL  TO  EXACT  AMBIENT  MAGNETIC  FIELD  MODEL 

RKM  ■ RE  A ALTF 
COLAT  s PIOV2  • ALATF 
ST  s SIN(COLAT) 

CT  m C08(C0LAT) 

SPH  m SlN(ALaNF) 

CPH  ■ C0S(AL0NF1 

calculate  magnetic  field  components  AT  POINT  P 
FROM  STASSINOPOULOS  M(10EL»  MODEL  5 (lORF  10/66), 


NEwmaG,S97 

NEmMAG,S96 

NEMMA6.S99 

NEmMaG.OOO 

NEMMAC.ttOl 

NCHMAG,<I02 

NEwMaG,«05 

nEhMAG.ROR 

nenmag.oos 

NEwMAC,-06 

NEnMA6,«07 

NEMMAG.aoS 

NEmMaG,«09 

NEmMAC.AIO 

NEHMAC.atl 

NENMAG,ttl2 

NEMMAG.ttlS 

NEMMAG.aiR 

NEHMaG.AIS 

NEMMA6,A16 

nemmag.oit 

NEnMaC.RIS 

NEHMA6.019 

NCmMaG,<I20 

NEmMAG.OZI 

NEmMa6,R22 

NEnma6,«2S 

NEnMaG.RZR 

NENMa6,02S 

NEmma6,«26 

NEMMA6,<t27 

NEmma6,«28 

NENMA6,«29 

NEMHAC.aso 

NCNMA6,ttSl 

NEmma6,432 

NEnMaG,43S 

NEMMAG,«3a 

NEMMA6,a35 

NEHMAG,a36 

NEmMa6,437 

NCmmaG,438 

NEnMaG,439 

NEmmaG,440 

NEkma6,441 

NCNMAG,442 

NEmMAG,443 

NEnMaG,444 

NEmMa6,445 

NE«Ma6,446 

NEMMA6.447 

NCmMaG,44S 

NEmmaC,449 

NEMMA6,4S0 

NC«maG,451 

NEmMA6,4S2 
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MAGFIT  (Cont’d) 


C*IL  nNEMG5(TM,RKM,8T,CT,SPM,CPH,fl«,BTHFTA,BPMl,B)  NEWM*G,«55 

NCwMAG.asfl 

calculate  magnetic  OIPOlE  momENT(GAUS8»KMJ)  from  the  field  CHMPOMENT8NCKHAtf,«55 

PR  (RADIAL),  BTHEtA  (POsITIvE  SOUTH),  ANO  BPHl  (PU8IT1VE  EAST).  NENM*6.«8fr 

NENMAG.aST 

R?SO  s BTHETAaBTHETA  ^ BPHI*bPhI  NENM*G,aS8 

HUO  * RKM**T*0.S*S0'>T(BP*8R  ♦ «t*R28a)  NENMA6«a59 

NEwMAG.ttdO 

calculate  the  sine  and  C08InE  OF  THE  ANGLE  ALPHA  AT  P WHICH  IS  THE  NEwHAG.AAI 
angle  between  the  magnetic  OIFOLC  HQHENT  and  geographic  north  pole  NEhMAG.RA? 

AND  DETERMINE  ThE  PROPER  SIGNS  NEmMaG.«6S 

NEWHAG.ttAtt 

ALPARG  « bphi/btheta  NEwMAG.OAS 

alpha  ■ ATAN(AlPARG)  NEwHAG.REG 

CnSALP  « AaS(CUS(ALPHA))  NCwMAG.aGT 

SINALP  c ABS(SIN(ALPha))  newhag.ras 

1F(RTHETA  .LT.  0.)COSALR  « •COSALP  NEwMAG.AfrB 

IF(BPM1  ,LT,  0.)SINALP  * aSTNALP  NEwHAG.ATO 

NEnHAG  ^71 

calculate  the  AnGlE  chi  measured  at  earth  center  which  NEwMAg|«72 

IS  the  angle  between  the  magnetic  dipole  moment  and  the  specified  NEwMAG.ATS 

POINT  p NEWMAG.ATR 

NEHMAG.ars 

CHIArG  m 2,*saRT(B28Q)/BR  NEwMAG,R76 

CHI  • ATAN(C'iIaRG)  NENNA6.tt77 

IFICHI  ,LT,  0.)CHI  * PI  ♦ CHI  NEWMAG,a7» 

CnsCHI  m C08(ChI)  NENMAG.a7« 

SINCHI  ■ StN(CHl)  NEWMAG.OBO 

NEnHAG.«81 

calculate  sine  and  COSINE  OF  THE  NEwMAG.aSP 

NORTH  LATITUDE  OF  THE  MAGNETIC  DIROLE  MOMENT  NEwHAG.ASJ 

NEwMAG.aBO 

SINLTO  • C08CHI*CT  ♦ 81NCHI*8T*C08ALP  nEwhaG,«S5 

CO8LT0  P 80RT(i,  • 8INLTo*SINLTO)  NEwMAG.OBS 

NEWMAC.UST 

calculate  the  east  longitude  of  the  magnetic  dipole  moment  NEwHaG.RBB 

NEmmaG.OBR 

STnOEL  ■ 8INCHI*sInAlP/Co3LT0  NEwHAG.RRO 

DEL  * ASIN(8INoEL)  NEwHAG.ARI 

COSDEL  c COSCHI  • SlNLTOwCT  NEwMA6,a92 

IFCCOSDEL  ,LT.  0,)0EL  « Sl6N(Pl,5INOEL)  • DEL  NEwMA6,093 

PHlO  s ALONE  • del  NEwMAG.aRO 

IF(PmI0  ,LT,  0.)  PHiO  * PHIO  ♦ TWOPI  NEwHA6,i!95 

NEwMA6,a96 

RETURN  NENMA6.497 

END  NEmMA6,R9S 
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SUBRUIITINE  nNf HCS(TM,RkM,9T,CT,SPH,CPH,BM,BTHETA,BPHI,B) 


* « • * 


THIS  RUUTTnF  calculates  the  macntic  ftelo  vector  at  a specified 

PnlNT  USING  mOPEL  % OF  STASSINOPOULOS  AND  nEAO  (nSSPC  72«12),  THE 
routine  is  a hopjFIEO  version  OF  ONEHAG  for  the  international 
GEUHAGNETIC  reference  FlEtO  (ICRF  IRbS.O). 


INPUT  parameters 


time  in  years  for  desired  field 
GEOCENTRIC  distance  OF  POINT  (MM) 

sine  of  geocentric  colatitude  of  point 

COSINf  OF  geocentric  colatituoe  of  point 

sine  of  geocentric  longitude  of  point  (POSITIVE  EAST) 

COSINE  OF  geocentric  LONGITUDE  OF  POINT  (POSITIVE  EAST 


OUTPUT  parameters 

BR  • radial  field  COMPONENT  (GAUSS) 

BTMETA  • POSITIVE  SOUTH  FIELD  COMPONENT  (6AU8S) 
BPHI  . POSITIVE  EAST  field  COMPONENT  (GAUSS) 

? • total  field  Magnitude  (gauss) 


dimension  LG(R«R)«LCT(9(9)«G(9,R),GG(9,9),GGT(9>9), 
1 8HMIT(R,R) 

DIMENSION  C0n8T(9»9),FN(R)»FM(R) 
dimension  P(9><>)«DP(9,9),SP(9)*CP(9) 


equivalence  (L6(ttI).66Cl*n)»(LCT(l,l)*&6T(l*l)) 


data  LO/1, 
A •203N»S05 
B ?«2,«176, 
C .R.lB.iaS 
D 13,-S,«ST 
data  LGT/i 
E -nef-tb? 
F ?9,.<I2.«2 
C «2««2r3t« 
DATA  SHMIT 
data  P(),l 


11 , J)/0,/»TMaLO/0,/»TiERU/l9bb,/»NH 

),CP(t).OP(l,I),SP(l)  f 2*1, ,2*0,  / 


,/»NhAX/9/ 


IF(8HMIT(1,1),E0,-I,)  go  to  8 

initialize  • ONCE  only,  FIRST  TIME  SUBROUTINE  IS  CALLED 


SHMlT(t,l)«*l, 

DO  te  N«l,9 
FN(N)sn 
DO  18  Mtl,9 
FM(M)sM«1 

IS  CONST(n,M)  c FlOAT((N,2)«*2«(M,1)**2)  / ( (2«N«3)*(?.*N»5) ) 
DO  ? N«?,9 

8HMlT(N,n  s (2*N.3)  • SHMlT(N.i,t}  f (N,l) 


JJs2 

DO  2 Ms2,N 


NEMMAG,a99 

NEmmaC.500 

NEmMaG.SOI 

NCmMaC.502 

NEmMaG.SOI 

NEMMAC.SOtt 

nehmag.sos 

NEmMaG,S06 

NEhmaG.SOT 

NEnMaC.SOB 

NEmmaC.509 

nEmmaC.SIO 

nenmag.sii 

NEmMaG.SIZ 

NEnmaC,S13 

)NEHHA6,Sia 

NEwma6,S1S 

NEmmaG.SU 

NEmmaG.SIT 

NEHMA6.5IB 

NCmHAG.519 

NEnMaC,S20 

NEnHa6,521 

NEHMAC4522 

NEmma6,S23 

NEnmaG,52« 

MEmma6,S2S 

NCmMaG,S2S 

NEmMaC,527 

NE*ma6,S2« 

NEmma6,S29 

nEnmaC,S3C 

»NEmma6,S31 

NEhmaC,S32 

NEmma6,S33 

NEMMA6,53a 

NEmHa6,S3S 

NEiiMAG,$3b 

nC»«maG,537 

NE‘iM46,S36 

NC*HAG,559 

nEhmaG,5«0 

NEnMAC.Sai 

nEnMaG,542 

nEmmaG,5«3 

NEMMAG.SttO 

NEhmaG,S45 

NEMMAG.Sttb 

NENMA6.Slt7 

NEmmaC,5«B 

NEMMAG.SR9 

NEmMAG.SSO 

NEmMaG.SSI 

NEmmaG.SS? 

NEMMAG.5S3 

NEMMAG.SSa 


i 


0NEMG5  ^Cont’d) 


& 

71 

SHHtTCN.Ml  « 8HMlT(fg,M.l)  * SORT  (ELOAT  ( (N>HA  1 ) * J J ) 

1 

107 

8HMIT(M.1,N)»8hMIT(N,M) 

g 

117 

2 

JJ  « 1 

'S 

12« 

71  ■ LCfl.n 

i 

125 

F2  « LGT(l,n 

127 

DO  7 Ns1»NMAX 

130 

DO  7 MalfNMAX 

140 

C6(N,H)  a Lr.(6,H)*8l*HIT(N,M)/ri 

1«S 

7 

GGT(N,M)  a LGT(N,M)*SMMtT(N»M)/F2 

if 

155 

8 

I7(TM,E0,TH0L0)  CO  TO  11 

\i 

% 

157 

TMOLOaTM 

160 

T a TM  . TZERO 

1 

161 

DO  10  Nal.NMAX 

163 

DO  10  HatfNMAX 

1 

173 

to 

6(N,M)  a GG(N,M)  * T*G6T(N,M) 

1 

g 

c 

r 

• *< 

k**  CAtCULATlON  USUAcLY  BEGINS  HERE 

1 

20« 

V 

11 

SP(2)aS2H 

n 

205 

C2(2)aCPM 

206 

00  12  Ma3,NMAX 

1 

215 

8P(M)a8P(2)*CP(M.nACP(2)*8P(H.l) 

220 

12 

CP<H)aCP(2)*CP(H«n«8P(?)*8P(H«l) 

1 

227 

A0Ra6371.2/RKM 

230 

ARaAOR**2 

1 

231 

BRaOyO 

232 

BTao.O 

1 

233 

BPa0,0 

/• 

23« 

on  21  Na2,NMAX 

2«2 

P(N«l,N)a0, 

1 

243 

21 

OP(N,l,N)ao, 

245 

00  17  NaZ^NHAX 

1 ■ 

252 

ARaAURaAR 

s 

253 

on  17  Hat,N 

i 

255 

ir(H,fO,N)  GO  TO  13 

1 

256 

I7(»I,E8,2)  GO  TO  19 

1 

262 

P(N,HlaCr*P(N<it,M)«C0N8T(N,M)ftP(N.2,M) 

265 

DP(N,MjaCT*0P(N«l,M),8T*P(N,j,M).CONST(N,M)*0P(N»2 

272 

CO  TO  14 

1 

301 

19 

P(N,M)aCT 

302 

OP{N,M)a«ST 

* 

300 

CO  TO  14 

se 

fS. 

307 

13 

P(N,N)a8T*P(N*l,N«l) 

311 

DP(N,N)aST«nP(N*l  >N«n«CT«PfN«l«N*l) 

i ■ 

517 

14 

PAR  a P(N,H)*AR 

323 

IP(H,E0,1)  GO  TO  15 

1 

332 

TFHPaG(N,M)«cP(M)tC(M«itN)*8P(H) 

t' 

335 

BPaHP-(C(N,M)*8P(H)«G(M*l»Nl«CP(H>)*7H(M)*PAR 

p- 

342 

CO  TO  16 

1 

346 

15 

TEhP  « C(N,M) 

352 

16 

BRaRR«TFilP«FN(N)*PAR 

557 

17 

BTaBTATEMP*nP(N,M)*AR 

571 

BPHI  a BP/8T/l00000t 

373 

BP  a BR/lOOnoO, 

1- 

370 

BTMETA  a Bt/100000, 

$ 

376 

B a SORT(BB*bP  ♦ BTHETA«8THETA  ♦ BPMUbPmI) 

404 

RETURN 

1 

oos 

END 
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nemhac.sss 

NEMM*G,5b6 
NExHAG.SS? 
NExMiG.SSe 
NEnMAG,S59 
NEmMaC.560 
NEnHa6,S61 
NEnma6,562 
NEmMaG.563 
NEmMAG«S6<I 
NENHAC.SbS 
NEkMaG.566 
NEmMa6,567 
NEn*4AC,S68 
NEhmaC.569 
NEMMA6.570 
NEnmaC,S71 
NENNA6.S72 
NEmMa6,S7} 
NEmmaG.57<i 
NENMA6.S75 
NEnha6,S7* 
nEi»maC,577 
NEnMa6.S78 
NEi«haG,579 
NENMAGtSaO 
NEnHaC.SSI 
NEMMAG.SaZ 
NEnHAK.SaS 
NExtiAG.Saa 
NEnHA(^,S6S 
NEMMA6.S86 
NEnmaC,S87 
NEKiMAG.saa 
NEMMAG.saa 
NEmmaC,S90 
»«EM'4AG,59t 
NEmMA6,S92 
NEi«Ma6.59S 
NEnMa6,S9« 
NEnHAG,S9S 
NENMA6.59b 
NEmMA6,S97 
NEf>MAG,S96 
NEMMA6.599 
NEi«MA6,bOO 
KE«»XAC,bOS 
NEk>4AG,602 
NEMMAG.bOS 
NEiiHAG,60« 
NEnMAG.AOS 
NEaHAG.bOb 
NE "MAG. 607 
NEhmaG,600 

NEWMA6.609 

NE»maG.610 

NENMAG.bll 

NEmMa6,612 
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Table  A8.  Summary  of  Running  Time  Experience  for 
Ambient  Magnetic  Field  Module 
on  a CDC  7600  Computer. 


Timing  runs  have  been  made  for  the  various  subroutines  in  the  ambient 
magnetic  field  model,  with  the  following  results  obtained  on  the 
Berkeley  CDC  7600  computer: 

MAGFIT  (includes  call  to  ONEMG5)  0. 30  msec^  or  0. 64  msec^ 

BFIELD  0. 055  msec 

CONJUG  0. 067  msec^ 

ONEMG5  0.  21  msec^  or  0.  56  msec^ 


For  a 6-page  Fortran  version  containing  no  DO-loops. 

“For  a 2-page  Fortran  version  containing  DO-loops. 

^This  number  should  be  contrasted  with  a value  of  26. 7 msec  required 
if  ONEMG5  (i.  e. , the  multipole  field)  were  used  instead  of  BFIELD 
(i.  e. , the  dipole  held)  in  tracing  the  field  line  to  the  conjugate  region. 
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